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Abstract 

The second moment of the heavy quark vector correlator at 0(a? s ) is presented. The 
implications of this result on recent determinations of the charm and bottom quark 
mass are discussed. 
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1 Introduction 



Correlators of quark currents are of prime interest for several phenomeno- 
logical applications. Their low-energy expansions, in particular, allow for the 
precise determination of charm and bottom quark masses via QCD sum rules 
[!T2ll3l4"ll5] . For this reason, heavy quark correlators have been frequently in- 
vestigated in the framework of perturbation theory. 

Up to 0{a 2 s ), analytic expansions to great depth are known for the low energy 
region. The three-loop QCD corrections to the correlator of two vector currents 
were first calculated in [B] . In [7] up to seven terms in the low energy expansion 
were obtained. This calculation also included further currents, namely the 
scalar, pseudo-scalar, and axial-vector current. Recently the calculation at 
three-loops has been extended to moments up to n = 30 for all four currents 
[81191. 



The moments of the vector correlator can then be used to extract the value of 
the masses of the charm and bottom quark from e + e~ data in the threshold 
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region using the i?-ratio, since they are related via a dispersion relation. A 
brief outline of this method is given in Section [21 which was first applied at 
three loops in [3]. 

At three loops a significant, sometimes dominant part of the error arises from 
the theoretical uncertainty due to higher orders, often estimated by the renor- 
malization scale dependence. Therefore the calculation had to be taken to the 
four-loop level [TOlfTT] to reach a precision comparable to or below the ex- 
perimental data. The contributions from double-fermionic loop insertions of 
heavy and/or light quarks are known explicitely up to 30 terms in the low 
energy expansion [T2]. The contributions due to light quark loop insertions 
of 0(a™n™ _1 ) are known to all orders in a s [T3]. Recently the lower moments 
were also calculated for the remaining three currents in p3]. 

In [1] the first moment of the vector correlator was used to extract the masses 
of the charm and bottom quarks. Since all but constant terms are known from 
renormalization group arguments, the analysis was done for up to the fourth 
moment, employing a conservative error estimate for the missing constant 
terms. 

In this paper we present the calculation of the second moment of the vec- 
tor correlator and discuss its impact on the determination of the charm and 
bottom quark masses. 

The outline of this paper is as follows: In Section [2] we set the framework and 
notations used throughout the paper. In Section [3] we explain the details of the 
calculation, present the result for the second physical moment and discuss its 
impact on the the quark mass determination. A brief summary and conclusions 
are given in Section HI 



2 Notation 



The correlator IP^g) of two vector currents is defined as 



with the current j**(x) = \I/(x)7 M \l/(x) being composed of the heavy quark 
fields ty(x). The function U liu (q) is conveniently written in the form 




(1) 



W u (q) = {-q 2 g» v + q^q v )Ii{q 2 ) . 



(2) 
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It can be related to the ratio R(s) = cr(e + e — > hadrons)/cr(e + e — > 
with the help of the dispersion relation 



127T Z jo s (s — q z > 



where the normalization 11(0) = has been adopted. 

To extract the quark masses the experimental data on the right hand side of 
(jnj) has to be compared with the theoretical evaluation of n(g 2 ) on the left 
hand side. This is best be done by comparing the corresponding Taylor series 
in q 2 . The n-th derivatives with respect to q 2 at q 2 = define the experimental 
moments 

M n P = I , (4) 



s 



which can be compared with the theoretical moments 

,9 / 1 



•^^iti^i c - (5) 



The latter are related to the Taylor coefficients C n of the vacuum polarization 
function 

fi(9 2 ) = ^IE^" (6) 

1U/I n>0 

with z = q 2 /(4m 2 ). Symbols carrying a bar indicate that the renormalization 
has been performed in the MS scheme. The coefficients C n can be expanded 
in a power series in ^ 

Cn = C$» + ^ + (^) 2 + (^) 3 Cg) + ■■. 

The four-loop contribution Cg> can be decomposed according to the number 
of quark loops and colour structures as follows: 

= CFTpn 2 C^ n + C F Tpn 2 h C^ n + C F T F nin h Ci^ n 

+ C F T F rii (CaC^a^ + CpC^n) + Cl (7) 

+ C F T F n h (C A C^ A n + CpC^l ^j + j^d abc d abc Csl . 

C contains the purely bosonic contributions, where we set the number 
of colours Nc = 3 for simplicity, while C Sn denotes the contribution from 

N 2 —1 

singlet diagrams. C F = 2 c Nc and C A = Nc are the Casimir operators of the 
fundamental and adjoint representation of the SU(N C ) group, respectively. 
Tp = | is the index of the fundamental representation. d abc is the symmetric 
structure constant, n\ and = 1 denote the number of light and heavy 
quarks, respectively. 
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3 Calculation and Results 



The diagrams have been generated using QGRAF [T5]- Expanding them in q 2 
results in four-loop tadpole integrals. Using EXP [16] they are mapped to six 
topologies with the maximum of nine lines. The main difficulty of the cal- 
culation lies in the reduction of the vast amount of integrals to the small 
set of 13 master integrals. This is done using Integration-By-Parts identities 
[2D] together with the Laporta algorithm [2T] which is efficiently implemented 
in the multi-threaded C++ program CRUSHER [17]. CRUSHER uses GiNAC [H] 
for simple algebraic manipulations and Fermat [19] for the simplification of 
complicated ratios of polynomials. A supplementary technique to perform the 
reduction to master integrals is based on the idea that self energy subgraphs 
of the integral can be reduced independently in order to effectively reduce the 
number of loops of the diagram. This can be useful because these integrals 
have up to two more propagator powers than integrals without an internal 
self energy and are therefore more cumbersome for traditional Laporta algo- 
rithm. In combination with Groebner Bases and the Mathematica package 
FIRE [22,23,21] it is also possible to calculate integrals without internal self 
energies. A more detailed description of the calculation techniques will be pub- 
lished soon [25J. In total the reduction of 1.8 million integrals was needed in 
order to perform the calculation, which is done using FORM [26J in combina- 
tion with the MATAD [27J setup. The necessary master integrals have been 
calculated in |28ll29ll3Uf31|l32[33|l34j . We confirm the results for the zeroth and 
first moment given in [8"|10|lll] . 

Inserting the master integrals and performing the renormalization of the strong 
coupling constant and the mass in the MS scheme leads to the following result 
for the second moment at fi 2 = m 2 as defined in Eq. (J7J: 



- (3) 64985074258811347 2900811008 

* " n_+ 353072079360000 3648645 ° 5 

1662518706713 , , 4 , , 2 , 362601376 , = 

(24a 4 + log 4 2 - 6C 2 og 2 2) + bg 5 2 

21016195200 v 6 s s ; 54729675 s 

725202752 , , 1684950406 , 

C2 log 3 2 U log 2 

10945935 s & 3648645 s 6 



725202752 , , 1684950406 , 

C2 log 3 2 U log 2 

10945935 s 6 3648645 s 6 

112680551036302633 26401638588211 164928917 

H 47076277248000 28021593600 270270 ' 

- (3) 5881974201847 97011619 , , 4 , , 2 N 

CSL = + + (24a 4 + log 4 2 - 6C2 log 2 2) 

5 ' 2 8369115955200 696729600 v B s & ' 

796232393699 745372259 

+ 371960709120^ ~ 185794560 ^ ' 

( 3) 20427854209619 31595849, 1 4n 2 . 

l°> ( 24a 4 + log 4 2 - 6C2 log 2) 

hNA,2 5649153269760 11612160 v & s s ; 

29638030087837 968787977 362 

~ 697426329600 ^ 3 + 15482880 ^ + ~63^ 5 ' 
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-O) 22559166733 520999 , , 4 , , 

C »« = - 10706100000 - 4354500 ( 24 " 4 + '° S 2 " 66 '° g 2 > 
309132631 167529079 

~~ 12902400 ^ + 5806080 ^ 4 ' 

- (3 ) 37320009196157 130387543 , , 4 > , 2 , 

Cri „ = (24o 4 + log 4 2 - 6C2 bg 2 2) 

ftA ' 2 271593907200 2177280 v 6 s 6 ; 

_ 5811074101069 2218910663 

6706022400 + 1451520 ' 
=0) 357543003871 520999 . 4o „. . 2 , 

C/7 9 = H 1 (24a 4 + log 4 2 - 6(2 log 2 2) 

lA < 2 11757312000 2177280 v ; 
36896356307 598455689 

~ 174182400 ^ 3 + 2903040 ^ 4 ' 
-( 3 ) 95040709 2029 , , 4 A , 2 N 

C '« = + 02765064 - 41472 ^ + >° g 2 " 66 '° S 2 > 
12159109 99421 

~ 4644864 ^ 3 + 55296 ' 
- (3) _ 1842464707 2744471 

hh ' 2 ~ + 646652160 ~ 1064448 ^ 3 ' 
- (3) _ 15441973 _ 32 

"> 2 19136250 45 



where Riemann's zeta function ( n and the polylogarithm Li„(l/2) are defined 
by 

OO 1 OO -y 

Cn = E^ and a n = Li n (l/2) = £-— . (8) 
fc=i K k=i z K 

For completeness we also give the results for the singlet contribution to the 
zeroth and first moment: 

- (3) 2411 6779 , 2189 , 5 , 73 / , An „ A , 2n 

Oi l = Cs + C4 Cs ( 24a 4 + log 4 2 - 6( 2 log 2 2 

5 '° 20160 4480 s 768 S 48 s 576 V 5 Si 5 



(9) 

w 3l 664837 2017831 , 175 > 739 / , 4 > , 2 x 

CL 3 = Cs + C4 ( 24a 4 + bg 4 2 - 6C2 bg 2 2) . 

5,1 2566080 855360 S 48 S 4320 V 4 ^ & Si s y 



(10) 

Numerically at fi 2 = m 2 one finds Cf ] \ ni=z = -3.49373 + 0.155877 and 
C 2 |n ; =4 = —2.64381 + 0.155877. The second term in each of these equations 
corresponds to the singlet contribution. 

Extracting the charm and bottom quark mass from the second moment using 
the input data given in [4] with the new value of leads to a shift of —3 MeV 
for m c and —2 MeV for mb and yields 

m c (3 GeV) = 0.976(16) GeV and m 6 (10 GeV) = 3.607(19) GeV . (11) 

This can be converted to the values at m c and m b , m c {m c ) = 1.277(16) GeV 
and rrib{mb) = 4.162(19) GeV, respectively. 

The final results for the quark masses given in [4j are rrib(mb) = 4.164(25) GeV 
and m c (m c ) = 1.286(13) GeV, respectively. In case of m c the first moment was 
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used at 0(a^) accuracy. For mb the second moment, which was known only 
up to 0(a 2 s ) at that time, was chosen. In the latter case the logarithms at 
0{a1) calculated by means of renormalization group methods were included 
and the error estimate was based on the missing constant term. Although this 
estimate was based on plausible arguments only a real calculation could prove 
its validity. Removing the 6 MeV error, which arises from the estimated term 
in case of the h quark, the total error of mb is reduced by ~ 25 %. In order 
«s the perturbative error is practically negligible and the remaining 19 MeV 
error arises from the experimental uncertainty and from the value of a s . At 
present this is the most precise determination of the bottom quark mass. 

As already discussed in [I], different moments weight the experimental re- 
sults from larger and smaller s values differently. Therefore it is important 
to compare the obtained quark masses from several moments to test the self- 
consistency of the method and the stability of the results. Because of sparse 
and poor experimental data in the continuum region above 4.8 GeV (for m c ) 
and 11.2 GeV (for m&), the data for R(s) were replaced by perturbative QCD 
in the analysis. This region can be suppressed by using higher moments, which 
is especially important in the case of m& where the first moment, which was 
already under full theoretical control at order in [I], receives a large con- 
tribution from the region above 11.2 GeV. The situation is significantly better 
for the second moment, which is now also fully under control from the theory 
side. For the determination of m c the first and the second moment are of equal 
reliability and the consistency between the two results for m c (3 GeV), namely 
0.986(13) GeV and 0.976(16) GeV, is remarkable. On the other hand for higher 
moments non-perturbative effects increase (especially for m c ) leading to larger 
theoretical uncertainties. For these reasons we think that for mb the second 
or maybe third moment are best suited for the mass determination, while for 
m c the first and second moment are preferred. 

Apart from the application discussed above, the higher moments evaluated 
above have been used recently for quark mass determinations from lattice 
simulations [5] and for the reconstruction of the full q 2 dependence of the 
vacuum polarization at C(af) [35J. 



4 Summary and Conclusion 

We have presented the second physical moment in the low energy expansion 
of the heavy quark vector correlator at four-loop order, including the singlet 
contribution. Although this contribution only causes a rather small shift in 
the quark masses obtained from the second moment the error is reduced sig- 
nificantly. The values remain in good agreement with those extracted using 
the first moment. 
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